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In a recent article we have studied the peculiar features of the Berry and Aharonov-Anandan's 
geometric phases for isolated spins S > 1. We have assumed that they are submitted to a dipole and 
quadrupole coupling to external E and B fields with the mild restriction E ■ B = 0. This implies 
discrete symmetries leading to remarkable simplifications of the geometry and algebra involved. The 
aim of the present work is to describe realistic proposals, within the realm of Atomic Physics, for 
the verification of some of our most significant theoretical predictions. There are several challenges 
to be overcome. For alkali atoms, most commonly used in atomic interferometers, the only practical 
way to generate quadrupole coupling, with a strength comparable to the dipole one, is the ac Stark 
effect induced by a nearly resonant light beam. One has then, to face the instability of the "dressed" 
atom hyperfine (hf) level, candidate for our isolated spin. One deleterious effect is the apparition of 
an imaginary part in the quadrupole to dipole coupling strength ratio, A. Fortunately we have found 
a simple way to get rid of ^(A) by an appropriate detuning. We are left with an unstable isolated 
spin. This implies an upper bound to the quantum cycle duration Tc. In the case of the Berry's 
phase, Tc has a lower bound coming from the necessity of keeping the non-adiabatic corrections 
below a predefined level. We have found a compromise in the case of the F = 2, m = ®''Rb ground 
state hf level. This is our candidate for the measurement of the somewhat "exotic" Berry's phase 
acquired by the 5* = 2, m = state at the end of a quantum cycle involving a rotation of tt of the E 
field - in practice the linear polarization of the dressing beam - about the B field direction. We have 
found a way to implement in a Ramsey-type interferometric measurement the procedures aiming at 
a control of the non-adiabatic corrections, as described in details in our previous theoretical article. 
A numerical simulation of our experimental proposal shows that a 0.1% accurate determination of 
Berry's phase, free of non-adiabatic corrections, can be achieved. Measurements could be considered 
also for cold ^^Cr chromium atoms with 5 = 3, where values of A ~ 1 can be obtained with an 
instability smaller than in the *^Rb case, due to a more favourable spectroscopic structure. The 
F = l,m — 1 hf level of the *^Rb ground state offers the opportunity to extend the measurement 
of Aharonov-Anandan's phases beyond the case S = |. We construct, using "light shift", the 
Hamiltonian H^^{t) generating a closed circuit in the density matrix space which satisfies at any 
time the "parallel transport" condition, thus making the quantum cycle free from the adiabaticity 
condition. We also consider the case of half-integer spins [e.g. ^"^Hg, ^"^^Ba and ^"^^Ba), with their 
own specific features. We show how the difference of Berry's phases for states S — ^ and S — ^, 
with m — ^, can be exploited to achieve an holonomic maximum entanglement of three Qbits. 

PACS numbers: 03.65.Vf, 42.50.Hz, 03.75.Dg, 37.25. -|-k, 03.67.Bg 



I. INTRODUCTION 



In a separate paper [l| we have presented a theoretical 
study of the Berry's phases generated by cyclic evolution 
of isolated spins of arbitrary large values. It was assumed 
that the spins interact non-linearly with time-dependent 
external electromagnetic fields (possibly effective ones) 
via the superposition of a dipole and a quadrupole cou- 
pling. We have made the assumption that the two ef- 
fective fields are orthogonal, a mild restriction but with 
many advantages. It implies several discrete symmetries 
of the spin Hamiltonian which simplify considerably the 
geometry and the algebra. Our purpose here, is to sug- 
gest atomic physics experiments to observe the original 
features of the Berry's phases that we predicted, but are 



still not revealed. 

In the present work we shall be mainly concerned by 
the adiabatic quantum cycles within a given time inter- 
val. They are generated by a Hamiltonian depending on 
a set of parameters assumed to be a system of coordi- 
nates for a differential manifold. In this way, an adia- 
batic quantum cycle generates a mapping of a close cir- 
cuit drawn upon the parameters space onto a closed loop 
upon the density matrix space. The Berry's phases can 
then be viewed as the geometric phases associated with 
this particular class of quantum cycles. 

The quantum mechanics postulates imply that Berry's 
phases can be written as Bohm-Aharonov loop integrals. 
The associated Abelian gauge field is acting within the 
space formed by the external parameters of the Hamil- 
tonian governing the quantum adiabatic cycles. In the 
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case of a non-linear Hamiltonian, involving both dipole 
and quadrupole couplings with the B and E fields sat- 
isfying the condition E • B = 0, the parameter space 
becomes isomorphic to the two-dimension (2D) complex 
projective space, CP2. Quite remarkably, CP2 can be 
identified with a solution of the Einstein equations in the 
4D real Euclidian curved space [3, Is)] . This kind of solu- 
tions appears in Quantum Theory of Radiation under the 
name of "Gravitional Instantons" Q . This strongly con- 
trasts with the magnetic dipole case where the parameter 
space is the familiar 2D sphere, and the Berry's phase 
gauge field the vector potential of a magnetic monopole, 
written in spherical coordinates. When both the dipole 
and the quadrupole couplings are present, the Berry's 
gauge field has a more complex structure, even for quan- 
tum cycles lying upon CP2 subspaces isomorphic to the 
2D sphere [23|- The observation of these non-trivial ge- 
ometry features, predicted by Quantum Mechanics, de- 
serves, in our opinion, a precise experimental investiga- 
tion. We suggest interferometric measurements, involv- 
ing light-shifted *^Rb hyperfine sub-levels, in order to 
exhibit these somewhat exotic quantum effects. 

It is quite natural to ask: how can one get from ex- 
periment the geometric phase associated with a quan- 
tum cycle of the density matrix p, since, by definition, 
p is phase independent? The answer is to be found 
within the Superposition Principle of Quantum Mechan- 
ics: any linear combination of two quantum states I'i'i) 
and [^'2) relative to a given quantum system, 1^1/12) = 
cil^'i) + C2 1^1/2) is an accessible state for the system. In 
the present atomic physics context, such a construction 
will be achieved via the interaction of the system with 
specific classical radio-frequency fields, using so-called 
Ramsey pulses 5]. The density matrix associated with 

1*12) is P12 = |*12)(^'l2| = IciPpi -f- |c2pp2 + Api2. 

The crossed contribution: Api2 = ci C2|5'i)(5'2| + h.c. 
contains all the information needed to obtain the differ- 
ence of the geometric phases acquired by the states 
and [^'2) during an adiabatic quantum cycle. We will 
discuss experimental schemes, where the geometric phase 
acquired by one state of the superposition will be known 
a priori to be zero, and one measures then directly the 
phase acquired by the second state modulo 27r. 

It is not possible to summarize here all the work stim- 
ulated by the original contributions of M. Berry, B. Si- 
mon, Y. Aharonov, D. Bohm and J. Anandan [a-Q- We 
refer the reader to review papers [III Ull) pedagogical 
T^- 1^ , several recent spin-off in quantum 
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ature only few papers dealing with the Berrry's phase 
for a spin submitted to a time-varying quadratic inter- 
action 24, 25j. The authors deal with the nuclear 
quadrupole resonance (NQR) spectra in a magnetic reso- 
nance experiment involving a rotating sample. However, 
they assume the absence of any magnetic interaction and 
this leads to level-degeneracy. Thus, the problem is gen- 
eralized to the adiabatic transport of degenerate states. 



In such a situation, the geometric phase is replaced by 
a unitary matrix given by the Wilson loop integral of a 
SU{n) non-abelian gauge potential where n is the dimen- 
sion of the eigenspace associated with a given degenerate 
quantum level [22i . i23l] . This makes an important differ- 
ence with respect to the conditions considered here, as 
well as in our theoretical work where level degeneracy 
was required to be absent. 

Throughout this paper we shall deal with a set of 
quadratic spin Hamiltonians 



HiBit), E(t)) = 75 S • B{t) + 7Q (S • m? 



(1) 



containing a Zeeman shift produced by a magnetic field 
B(t) (real or effective) and a quadratic Stark shift pro- 
duced by an electric field E(i) (mainly effective). The 
non-linear spin coupling is responsible for new physical 
features becoming apparent for S > I. 

A crucial step, in the confrontation of the theoretical 
results with experiments, was to realize that, rather than 
the standard dc-Stark effect, the ac-Stark shifts, induced 
by a nearly resonant linearly polarized laser beam, were 
the proper tools to generate an effective E field. Samples 
of cold alkali atoms or trapped alkali-like ions appear 
then as good candidates for the observation of the Berry's 
phases induced by non-linear coupling. Our spin system 
S is identified with the total angular momentum F acting 
upon a given hyperfine (hf) sub-level of an alkali atom 
ground state. 

Our method relies upon the second-order ac Stark 
shifts involving optical frequencies nearly resonant for 
the transitions 6*1/2, f — ^ Pi/2.F', with an appropriate 
choice of the detunings [2^. This would then allow to 
use the hyperfine ground state sub-levels of rubidium and 
cesium isotopes to simulate isolated spin systems with S 
having integer values between one and four. In addition, 
the laser frequency can be tuned for making one of the 
hyperfine sub-level insensitive to the laser field. It will 
then provide an absolute phase reference for the second 
hyperfine sub-level which is performing a quantum cycle. 

However, the method of ac-Stark shift to construct the 
quadratic coupling, presents one drawback. Under the ef- 
fect of irradiation by the light field, the hyperfine ground 
state sub-level used to simulate our spin system will ac- 
quire a finite decay rate. This puts an upper limit upon 
the time available for performing one quantum cycle and 
detecting the associated phase shift. There are atomic 
species for which this detrimental effect is much less se- 
vere. For instance the ^°^Hg mercury isotope, (and alkali- 
earth-like atoms with half-integer nuclear spin) have a 
spectroscopic structure more favourable for this purpose. 
It appears that ^^Cr chromium atoms with a spin 3 of 
purely electronic origin could provide suitable and espe- 
cially interesting candidates. 

The derivation of the Berry's phase requires the va- 
lidity of the adiabatic approximation. This means that 
to perform a valid measurement, one must exert a very 
tight control upon the non-adiabatic corrections, which 
are governed by the Hamiltonian-parameter velocities. 
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This problem is particularly crucial, here, because of the 
instability introduced by the ac Stark shift in our spin 
system. It is addressed in great detail in our previous 
theoretical work T] . We have proposed explicit solutions 
for approaching the adiabaticity criterion while ramping 
up the E field and applying the angular rotation speed of 
the periodic Euler angles. We have given definite proce- 
dures for reducing, eventually suppressing non-adiabatic 
corrections of various origins. 

As two examples of the new features of the Berry's 
phases generated by the quadratic Hamiltonian let us 
consider two particular set of adiabatic cycles within the 
time interval < t < T^. 

i) The B field is precessing around a fixed axis, while 
the E field (orthogonal to B) is lying within the rotating 
plane defined by B and the rotation axis. 

ii) The direction of the B field is fixed and the orthogo- 
nal E field is rotating around B by an angle a{t) with the 
boundary condition a{Tc) — a(0) = tt modulo tt. Per- 
forming an "adiabatic" increase of the E field starting 
from a null value, we have obtained in [l| a mathemat- 
ical form of the Berry's phase relative to cycle i) which 
looks superficially similar to that of a linear Hamiltonian. 
However, there is one significant difference : the contri- 
bution involving the cosine of the B field tilt angle is no 
longer proportional to the magnetic quantum number m 
whatever the spin value, but instead to the spin polar- 
ization along the B field. In addition, a new contribution 
to the Berry's phase is generated at the end of the cycle 
ii). It is given by the spin polarization times the rotation 
angular speed a{t) integrated over the < t < Tc inter- 
val. A stricking case is that of m = for a spin larger 
than 1, the usual Berry's phase generated by the linear 
Hamiltonian has a null value, while the one generated by 
the quadratic Hamiltonian of Eq. ([TJ is non-vanishing 
and increases with the magnitude of E up to a maximum 
growing with the value of the spin (especially when it is 
an even integer). Our aim, here, is to define the precise 
experimental procedure required for the observation of 
such new features of the Berry's phase, emerging from 
our previous work 

We have already underlined that the E, B field orthog- 
onality condition endows our problem with important 
symmetry properties and that it also makes the param- 
eter space isomorphic to CP^, but in addition, for the 
special case of spin one, this condition implies the isomor- 
phism of the density matrix with the parameter spaces. 
Thus, the Berry's phase generated by the quadratic spin 
Hamiltonian was found to be mathematically identical 
to the Aharonov-Anandan (AA) phase using an appro- 
priate parametrization of CP^. However, the physical 
contents are in general different, since, in contrast to the 
Berry's phase, the AA phase is not restricted to adiabatic 
quantum cycles. This case was previously discussed the- 
oretically by C. Bouchiat and G.W. Gibbons [l^ and 
C. Bouchiat [2^. Using cold ®''Rb atoms in the hf state 
_F = 1 experiencing a quadratic ac Stark shift, we propose 
here to perform an experimental comparison between the 



AA and Berry's phase in different adiabatic and non- 
adiabatic regimes, first, to verify, within the adiabatic ap- 
proximation, the identity between Berry's and AA phases 
for appropriate parametrization and second, to observe 
the difference in their behaviour when the adiabatic ap- 
proximation begins to fail. We show how a geometric 
phase equal to Berry's should be still obtained by per- 
forming a quantum cycle within the condition of "paral- 
lel transport" , Tr{p{t) H\\(t)) = 0, without any constraint 
upon the time derivatives of the physical observables. We 
shall discuss the realization of measurements of this kind 
with ®^Rb atoms in their hf state F — 1, showing inter- 
est for the fast accomplishment of elementary operations 
in quantum processing. For spins larger than one, the 
Berry's phase relative to the Hamiltonian considered in 
this work loses the identity with the A A phase, which in- 
volves now closed circuits drawn upon larger projective 
complex planes CP"^^ . 

Finally, as an additional support to the present exper- 
imental program, we would like to point out that it may 
have other significative physical spin-off. A first simple 
suggestion would be to stop midway the adiabatic cy- 
cle described in the present paper. For an appropriate 
choice of the Hamiltonian parameters, one obtains an 
example of "coherent spin-squeezed" states (30l - [3^ . A 
second example relies upon the non-trivial dependence 
upon of Berry's phases of the present paper, in con- 
trast with the dipole case where Berry's phase is pro- 
portional to the magnetic quantum number m whatever 
S. This property is the basic ingredient used to perform 
an "holonomic entanglement" of non-correlated spins 
1/2 (or Q bits) having a fixed number of spins "down". 
The corresponding vector state can be written as a lin- 
ear combination of the eigenstates of and Sz 

where S — J^iLi the total spin operator. Despite 
the fact that a given eigenvalue of may appear sev- 
eral times among them, the states 5'^ can always be 
chosen to have different symmetry properties upon the 
permutations of the TV spins, insuring their orthogonal- 
ity. If we assume that in formula (1) S stands for the 
total spin operator introduced above, the corresponding 
Hamiltonian H]\r{t) is clearly invariant upon all the per- 
mutations of the N spins. As a consequence, it acts upon 
the states ^'^ j^.j as if they were isolated spin systems. At 
the end of the Berry' cycle, they acquire a phase depend- 
ing upon their S values. This implies an entanglement of 
the initially non-correlated N spin states f^. In the two 
particular examples with = 3,4 and one spin down, 
i.e M = S — 1, the parameters of the cycle can be cho- 
sen in order to achieve a maximum entanglement. (The 
simplest non-trivial example, A^ = 3, is discussed in the 
Appendix). 

Sec. II is a summary of our theoretical work 1,] which 
presents the remarkable properties of Berry's phases ac- 
quired by an arbitrary spin, non-linearly coupled, at the 
end of a closed adiabatic quantum cycle. We also re- 
mind different procedures for keeping the non-adiabatic 
corrections below a predetermined level. In Sec. HI, the 
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expression of the quadratic spin Hamiltonian Eq. ([T]) 
is derived in terms of the experimental parameters for 
alkah atoms lighted up with a laser beam close to reso- 
nance. We present new possibilities offered by cold ^^Cr 
chromium atoms which are now available. Sec. IV de- 
scribes a Ramsey type interferometry measurement of 
Berry's phase, free of any non-adiabatic correction, ac- 
quired by the F = 2,m = hf substate of ^^Rb, at the 
end of a quantum cycle induced by an E field rotating 
around a fixed B field. Sec. V underlines the pecu- 
liar features of half-integer spin Berry's phases, also yet 
non-observed. They are illustrated by the case of spins 
three-half, where two types of measurements look possi- 
ble. One involves ^°^Hg atoms (or other alkali-earth-like 
odd isotopes) , the second a nuclear quadrupole resonance 
experiment on a uniaxial crystal placed in a magnetic 
field. Sec. VI deals with the particular case S=l. We give 
a method to achieve parallel transport on *^Rb F = 1 
atoms and perform an experimental comparison between 
Berry's and AA phases. Finally Sec. VII is a summary of 
our work and possible development. 



II. SUMMARY OF BERRY'S PHASE THEORY 
FOR ARBITRARY SPINS NON-LINEARLY 
COUPLED TO EXTERNAL FIELDS 

Throughout this paper, we make two important as- 
sumptions, first the two spin couplings, linear and 
quadratic, can be made of comparable magnitudes, sec- 
ond the effective fields B and E are orthogonal. 

A. Symmetry properties of the non-linear spin 
Hamiltonian. Physical implications 

In order to exhibit the geometric structure of 
the Hamiltonian i/(B,E), let us introduce the rota- 
tion R{t) mapping the coordinate trihedron (x,y,z) 
upon the trihedron defined by the directions of the 
B,E fields: R{t) = H{z,ip{t))H{y,e{t))H{z,a{t)), 
where TZ{u, x(^)) is standing for the rotation of an- 
gle x(t) around the unit vector u. A very conve- 
nient mathematical tool is the unitary transformation: 

U{R{t)) = exp-^SMt)exp-j-Sye{t)exp-j-SMt), 
which performs the rotation of the spin operator: 
WiRit))SU{R{t)) ^ R{t) ■ S. This rule allows us to 
write: i?(B,E) = 7sS?i U {R{t))n{X)W {R{t)) , with 
'H{\) = Szh~^ + XSlh^^. Taking aside the energy scale 
jsBh, the dimensionless parameter A = Ti^qE'^ / (^sB), 
combined with the Euler angles 0,Lp,a provides a set of 
dimensionless parameters for H(B, E), which could serve 
as coordinates for the CP^ projective space. 

The eigenstates of H(A), tp{X,m) are labeled with a 
"magnetic" number m by requiring that their analyti- 
cal continuation towards A = coincide with the angu- 
lar momentum eigenstates | S, m) . Similarly, the associ- 
ated eigenenergies £{\,m) satisfy the boundary condi- 



tion: £"(0, m) = m. 

Thanks to the constraint E • B = 0, 'H(A) has sev- 
eral discrete symmetries, leading to important physical 
consequences 

• a. The "m" parity (—1)'^"™, associated with a tt- 
rotation around z, is a good quantum number for 
■H(A). As a consequence, 'H(A) cannot mix angular 
momentum states l^*, m) with opposite "m" par- 
ities and it can be expressed, within the angular 
momentum basis, as the direct sum of two matri- 
ces acting respectively upon the states with even 
and odd m parity: H(A) — Heuen(A) HoddW- 

• b. 'H(A) obeys under the rotation TZ{x, tt) the trans- 
formation law: H{X) — A). Since the rota- 
tion TZ{x,'k) fiips the spin component along the z 
axis, one gets the symmetry relation: £{m,X) = 
— £(— m, —A). 

• c. The invariance of 'H(A) upon the rotation 
TZ{z, tt) implies that the quantum average of S rela- 
tive to the vector state ip{X, m) lies along the z axis: 
{4>{X,m)\S\'ilj{X,m)) = hp{m,X)z where p{m,X) 
will be referred as the polarization of the spin state. 
Remembering that the unit vector z is taken along 
the B field, one gets for the following quantum spin 
average: ( S ) = ?ip(m, A) B/B. 

• d. The last invariance is somewhat more subtle 
than the three others since it involves the antiuni- 
tary transformation associated with the product of 
the time reversal by the space reflection with re- 
spect to the xy plane. It implies that ip{X,m) can 
be represented by a real vector. As we shall see this 
result has non-trivial physical consequences. 

The eigenvalues of i?(B, E) together with the associated 
eigenvectors \I'(m, A) can be written as: E{m,B,E) = 
jsB £{m, A) ; ^{m, X,t) = U {R)tp{m, A). The spin quan- 
tum average relative to 5'(rn, A) is along the B direction: 
(S) = fip{m, X)'B/B. Using the Hellmann-Feynmann 
theorem, p{m. A) is obtained by taking the partial deriva- 
tive of the energy E{m,B,E) with respect to 7s-B: 
p{m,X) = £{m,X) — A^^^^^. Note that p{m,X) obeys 
under the reversal of A the same symmetry law as £{m. A) 
derived in • b. The curves representing the variations of 
£{m, X) and p{m, A) versus A for S = 2,3 and 4 are dis- 
played in Fig. 1 of our previous paper [l| . 

B. A sketchy derivation of Berry's phases for 
adiabatic cycles governed by ff(B,E) 

We have now all the necessary ingredients to sketch the 
derivation of the Berry's phase using the instantaneous 
eigenfunctions ^(m, A,i) given in the previous subsec- 
tion. Our starting point is the standard formula: 

/■^ d 
/3(m) = J dt{^im,X,t)\i—^{m,X,t))+(l}{m) 
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FIG. 1: (Color online) Reduced energies £{m, A) (left-hand) and gauge field components Ac — pirn, \) — m (right-hand) versus 
A for S — 2 and S = 3, for selected m-values relevant for the discussion presented in subsec. 2.C. Intersection of the energy 
curves with the vertical axis indicates the eigenvalues m of Sz for A = 0. For A > remarkable effects appear for S = 2, m = 
and for S — 3, m = —1. In both cases this occurs in a range of Stark to Zeeman coupling ratios where the levels are well 
separated. 



with (^(m) = arg(*(m, A(0),T)/*(m,0,A(0)), (2) 

where we have used the fact that A is a non 
periodic parameter: A(T) = A(0). Inserting 
\I'(m, A,i) — U{R)ip{ni, X), the time integral contri- 
bution can be rewritten as: 

7(m) = dtzhii^im, A)| (f + W{R)§-^U{R)) i>{m, A)). 

Since ■(/'(to, A) is a normaUzed real vector, the first term 
inside the parenthesis vanishes. In the second term, the 
operator i UW (R) -§^U (R) is a standard group theory 
object which can be expressed under the canonical form: 
iKU\R)§-^U{R) = uj{t) ■ S, where uj{t) is a real vector. 
Its three components are linear functions of the time 
derivatives of the Euler angles. Their explicit expressions 
can be found in references [l^. Using • c, one sees im- 
mediately that only the z- component Wz = cos{9)^p -\- a 



does contribute: "f{m) = dtp{m, X) {cos{9) ip + a). 
The calculation of (j){m) is performed in details in 
references It is greatly simplified by the particular 
features of the expansion of ^{m, A) over the angular 
momentum eigenstates resulting from properties • a. 
and • d, which lead to : 

|^(TO,A(0)) = Eim-2n\<sCmAm)\S,m - 2n). The 
sum runs upon states of the same m-parity and the coeffi- 
cient Cm,n are real numbers. One can factorize out from 
C/(i?(T))C/t(i?(0)) the operator giving the phase shift 
0(m): = exp-^5,(^(r)-(^(0) + a(T)-a(0)) 

Remembering the quantum cycle boundary conditions: 
ip{T) — (p{0) = 2n^ TT , a{T) — a(0) = Ua tt, one sees 
that the effect oi upon the terms of the expansion of 
'0(m, A(0)) is just to multiply them by the same phase 
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FIG. 2; (Color online) Plot of £{rn,\)/\ within the interval 0.1 < A < 2, allowing us to clarify the rearrangement of 25' + 1 
non-degenerate levels for A ^ 1, into S degenerate doublets with £{m,X)/X ~ fj,^ in the limit X ^ 1, fj, being an integer such 
that —5* < j-L < S. Starting from the even-odd (or odd-even) pair: {£{S, X)/ X, £{S — 1, A)/A}, one sees clearly that the pair 
converges, without crossing, towards the degenerate doublet with: £{S,X)/X ~ £{S — 1, A)/A ~ S^. The next lower pair will 
end as the degenerate doublet having the energy ~ A (S* — l)'^ and so on, until one reaches the isolated level with /i = —S. It 
has no other possibility than to converge to the non degenerate level with fj/' — 0. The fact that there is no level crossing for 
finite values of A follows from the very simple mathematical structure of the Hamiltonian 'H{X). Indeed, any modification of 
its symmetry properties, is excluded until one reaches the limit A — >■ oo. 



factor: exp-im{ip{T) ~ ip{0) + a{T) - a{0)). In this 
way one obtains (/)(m) — — m{ip{t) + a{t)) dt. 

We arrive at our final expression for the Berry's phase 
in terms of a loop integral along a closed circuit drawn 
upon the parameter space CP^: 



l3{S,m) 

A^{S,m,X) 
Aa{S,m,X) 



A^{S, m, A) dip + Aa{S, m, A) da, 



p{S, m. A) cos 9 
p{S, m. A) — m. 



m, 



(3) 



In the writing, we have stressed the two new physical fea- 
tures introduced by the quadrupole spin coupling: first, 
the Aharonov Bohm-like integral now involves a two- 
component Abelian gauge field {A^,Aa) instead of a 
single one, second, it now exhibits a strong dependence 
upon the value of = h'^S{S + 1), in contrast with the 
dipole case where the Berry's phase depends only upon 
m. The latter effect refiects the fact that the polariza- 
tion p{S, TO, A) has a simple linear relation with £(rn, A). 
These eigenenergies are given by the roots of two polyno- 
mials having respectively the degree 5-1-1 and S —1 and 
coefficients which are A monomials. As a consequence, 
for 5* > 4 the eigenenergies are given by transcendental 
functions of A. 



C. Physical implications of the S-dependence of 

Berry's phases generated by iif(B,E) 

For integer spin values, remarkable effects are pre- 
dicted when Berry's cycles have as initial state the an- 
gular momentum eigenstate jS', 0), corresponding to an 
initially vanishing polarization p{S, 0, 0) = 0. If no E 
field is applied during the cycle, the polarization keeps 
its null value all along the cycle and one recovers the well 
known result I3{S, 0, 0) = 0. However, if the E field in- 
tensity is ramping up to a maximum before a-rotation 
starts and returning to a null value when it stops, the 
inital state is mixed with states \S,m) with m ~ ±2n 
(n is an integer < S/2). For 5 > 1 a finite polar- 
ization p{S, 0, A) appears to third order in A given by 
p(0,A) = ^X^S{S + 2){S^~l){l + OiX^)). Finite values 
for the gauge field {A^,Aa) are then generated. 

In the following we concentrate upon the cycles where 
a and A are the sole varying parameters. As a con- 
sequence, only the gauge field component Aa is rele- 
vant. We have displayed in Figure 1 the variations with 
A of £{S,m,X) for 5 = 3 and 2 with to = 0,±1,±2, 
together with that of the gauge field ^^(S', to. A) for 
TO = 0, —1, —2. One sees clearly that the relevant energy 
levels are well separated. The magnitude of Aa{S;m, X) 
is remarkably large for S" = 2,to = 0, in strong contrast 
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with the case of a hnear spin Hamihonian. We note that 
Aa{S,0) is smaller for 5' = 3. For |A| < 1, this follows 
from the fact that the term cx has a large coefficient 
which changes its sign with the parity of S. For to = — 1, 
the largest value of AaiS,—!) occm'S for S = 3. We 
note that p(3, — 1,A) > when A > 0.7. This indicates 
that the state |3, — 1) is strongly mixed with the states 
|3, TO = 1, 3) within the range .7 < A < 2. 

To conclude this subsection, we emphasize that the 
holonomic entanglement method of > 2 non-correlated 
^-spins, described in reference relies upon the fact 
that Berry's phase has a strong dependence upon S. A 
different example is given in Appendix A. In a theoret- 
ical note we show how to exploit the difference between 
Berry's phases for S — ^ and S — ■^j with to = ^, to per- 
form an holonomic entanglement between three initially 
non-correlated ^-spins. 

D. A summary of the non adiabatic corrections to 
Berry's cycles 

To perform empirical determinations of the Berry's 
phase, one must tackle the problem of the non-adiabatic 
corrections. This question becomes crucial in the exper- 
imental situations where the coherence decay time in- 
duced by application of the Hamiltonian i?(B, E) puts 
constraints on the cycle duration. Shortening the quan- 
tum cycle risks to spoil the validity of the adiabatic ap- 
proximation. In the present context, we have found con- 
venient to study the quantum cycle in the rotating frame 
attached to the time-varying fields. The Coriolis effect 
generates an extra magnetic field AB which involves a 
linear combination of the Euler-angies time derivatives. 
The longitudinal component along the B field is the only 
one which survives when a is the sole time-dependent 
Euler angle. The rotating frame Hamiltonian H^^ is 
void of any geometry. As a consequence, the phase shift 
acquired at the end of the cycle (/)(B -I- AB) is purely 
dynamical. The laboratory Berry's phase contribution 
7 (to) is incorporated into the dynamical phase under 
the form of its first-order contribution with respect to 
rj = —AB///B = {cos6 (p + a)/{'jsB). The higher-order 
terms give all the non-adiabatic corrections associated 
with a when it is the only varying periodic parameter. 
We have also shown that the subset of these corrections, 
odd under a reversal of 77, cancel exactly for "magic" val- 
ues A = A*(ry), obtained from a polynomial fit performed 
on the numerical results valid for < 77 < 0.5 (see P 
Sec.V.A) 

A* (2,77) = 0.838213 - 0.0837823 77^ _ 0.0431478 77-* - 
0.0231887 77*^ - 0.0207986 77*^. (4) 

This cancellation is implemented in the experimental 
project described in Section IV. The 77-even corrections 
are eliminated by subtracting the phases measured for 
two "mirror" cycles (77 — >■ —77). The case of the non- 
adiabatic corrections induced by the transverse field 



AB±/B is somewhat more involved since it introduces 
a non-trivial geometry and, as a consequence, a Berry's 
phase contribution to be added to the one coming from 
the transverse dynamical phase; the results explicited in 
[l| are not used in the present context. 

For the adiabatic approximation to be satisfied, the 
primary condition is that, at each instant, the spin quan- 
tum state is an eigenstate of i?(B(i), E(i)). In the above 
discussion, we have made implicitly the two following as- 
sumptions: i) during an Euler angle cycle, X{t) has a 
slowly varying value of the order of unity, ii) the ini- 
tial state of the Euler cycle is obtained by an adiabatic 
ramping governed by ■H(A(i)) = Sz + X{t)Sl, starting 
from a null value of A in order to get its desired value 
for a well defined value of to. We have shown, by doing 
explicit calculations, that approaching this ideal "adia- 
batic ramping" is, in practice, a non-trivial task. The 
key-parameter which governs the time dependence of the 
quantum state is the time derivative of A. A linear in- 
crease of A(t) would be equivalent to a rf pulse with sharp 
edges, leading to large oscillating non-adiabatic correc- 
tions exhibited in ref. [l|. Among several methods [33j . 
a standard procedure to smooth them out is to use a 
Blackman pulse shape [l^]. In the present context, this 
condition is implemented by taking for X{t) the following 
time dependence: A(i) = Xof{t/T)/f{0), where f(s) is 
the Blackman function 

/(s) = 0.42 - 0.5 cos {2ns) + 0.08 cos (47rs) (5) 

and T the ramping time. The efficiency of the procedure 
has been illustrated before (see Fig. 5 in [l[). 

There is a third assumption implicit in the rotating 
frame analysis, namely that the adiabatic approximation 
is valid for the rotating frame Hamiltonian i?//. As will 
be shown by the theoretical analysis of the experimental 
project presented in section IV, the adiabatic approxima- 
tion works beautifully provided one uses also a Blackman- 
pulse shape for the angular speed, a{t). 

When looking at a plot of £{m,X)/X (Fig. 2), it ap- 
pears clearly that there are pairs of states associated 
with TO values differing by 1 which never cross but be- 
come quasi non-degenerate for large values of A. When 
one is interested in quantum cycles where only a and 
A are time-dependent this should not affect the results 
since the Hamiltonian _ff (B(i), E(i)) has no matrix el- 
ement connecting the Ato = ±1 states. Nevertheless, 
when two such levels happen to be close, the spin system 
becomes particularly sensitive to imperfections which al- 
ter the symmetry, as for instance a stray component of 
the magnetic field orthogonal to B. 

III. ATOMIC SIMULATION OF ISOLATED 
SPINS NON-LINEARLY COUPLED TO 
EXTERNAL FIELDS. 

Since, among cold atoms, alkali are the most frequently 
studied we first examine what kind of measurement looks 
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possible in their case. Such a choice is motivated by 
the fact these systems can be kept in a decoherence-free 
space for a relatively long time (up to 1 s), if cooled and 
trapped. But as is well-known, it is very difficult to ob- 
serve a static quadratic Stark effect on hf sub-levels of al- 
kali ground states, since their electric tensor polarizabil- 
ities are strongly suppressed. A tiny effect appears only 
as the result of a third-order perturbation 35, 36] , when 
the tensor part of the hf interaction is taken into account. 
The tensor polarizability of Rb 012 —2 mHz/(kV/cm)^ 
[131, implies that a field of w 300 kV/cm would be re- 
quired to generate a coupling strength of 100 Hz. This 
makes it unrealistic to use a rotating static E-field to 
observe the quadratic Berry's phase. 

On the other hand one may rely on light shifts "2^. 
Hereafter, we consider the case of a quadratic coupling 
induced by linearly polarized light fields which have the 
advantage, with respect to microwave- or rf-fields, of 
making the orientation of E relative to B both precisely 
ajustable and easily rotated. However, with this method 
one has to face the problem of instability of the dressed 
atomic state if one wants to obtain a quadrupole to dipole 
spin coupling ratio of order 1 . This is not straightforward 
for the ground state hyperfine levels of alkali atoms. In 
this case, the role of the spin operator S is played by the 
total angular momentum operator ?iF = ?i(s -f I), where 
Tis and fii are the electronic and nuclear spin operators. 
Actually, we are going to show that the practical real- 
ization of the Hamiltonian i/(B(t), E(t)) required to test 
the theoretical predictions of this work for A ~ 1 appears 
possible for *^Rb atoms. 

Now, in the rapidly expanding family of laser cooled 
and trapped atoms has appeared the daring alternative 
of chromium. The great progress achieved in the manip- 
ulation of this atom ha ve g iven rise to a series of beautiful 
experiments, t.g. |38| - |40| |. As a spin S* = 3 candidate, 
it has the advantage of being of pure electronic origin. 
This makes the instability of the dressed atom no longer 
a problem. The chromium option is discussed at the end 
of this section. 



A. Building quadrupole spin couplings in alkali 
atoms using the ac Stark effect. 

It is known that the application of a light beam close to 
resonance with one atomic excited state generates light 
shifts which can simulate the effect of an electric field 
or a magnetic field [26|. A fictitious B field arises if the 
beam is circularly polarized while a fictitious E-field is 
created by a linearly polarized beam. One might won- 
der therefore, whether, with a single beam of elliptically 
polarized light, it would be possible to generate both ficti- 
tious electric and magnetic fields satisfying the condition 
E-B = 0. We consider the case of alkali atoms irradiated 
by a laser beam nearly resonant with one hf component 
of the nSxji p — nPxji jr transition. The atom-laser cou- 
pling responsible for ac Stark shifts - or light shifts - of 



the ground state sublevel can be calculated to second 
order in the atom-radiation field interaction, in the ro- 
tating wave approximation (see [1^, If we make the 
simplifying assumption F = !F , discussed hereafter, we 
can use the proportionality between the atomic electric 
dipole and the angular momentum operators, resulting 
from the Wigner-Eckart theorem, D/eao = dgpF, where 
d is the Arris = matrix element of this electric dipole 
transition in atomic units and gp = 2(i^ — /)/(/ + l/2). 
In this case, the light shifts of the hf ground state Si/2^f 
can be represented very simply in terms of the effective 
Hamiltonian 

HisiF) = —gl X 

^ie*Ae •Ffc + i(F-e* F-e + F-eF-e*)^ , (6) 

where = d£ /2h represents the dipolar coupling of the 
atom with the laser field. The field magnitude S is related 
to the photon number density by e^E"^ ~ Nfiuj/V [il], 
[42! and the complex vector e of unit norm defines the 
polarization. For simplicity, we suppose the detuning 
A between the laser and the atomic transion frequency, 
i-^F.j^, to be small compared to the hf splittings, so that 
we can consider the contribution of the J-' — F line 
alone, A = cj — (jJf,t- In the next subsection, we shall 
consider the case where this condition no longer applies. 
It will turn out that the general form of His{F) is the 
same but with different coefficients. 

A general expression for the elliptical polarization e of 
a beam directed along z is: 

e ~ —j= (cos (5e e_|_ -I- sin (5e e_ ) 
v2 

where e± = -j=— exp (=F ru), (7) 

v2 

and the angle u defines the orientation of the ellipse in 
the (x, y) plane of polarization. 

For rr = the expression for Fli^ obtained after some 
angular momentum algebra is: 

HisiF) = -^gl X 

f „ 9 , . > 9s 1 — sin 25 e \ 
{ cos 2(5e F, + sin 26e F^ + {F{F + 1) - F^) ^ J 

(8) 

The first term, linear in F involves the circular polar- 
ization of the beam, characterized by its helicity ^ — 
Im{e* A e • fc} = cos2(5e, while the quadratic contri- 
bution cx F^ involves the linearly polarized intensity. 
Contrary to a static E-field, a laser field can provide a 
sufficiently strong quadrupolar coupling, for laser power 
and frequency adjustments lying within convenient limits 
(see Sec IV). However, because the ac- field is complex, 
while the dc field is real, His{F) actually differs from 
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the Hamiltonian H{B,E) of Eq.(IT]) by its last term in- 
volving , which cancels out only if the polarization is 
purely linear, i.e. S^, = 7r/4. Although this does not ren- 
der the light-shift Hamiltonian untractable, one looses 
the simple symmetry property of the Hamiltonian under 
the rotation TZ{x,tt) discussed in section H.B. For this 
reason, to keep our concrete discussion of the theoretical 
implications as simple as possible, we prefer to choose an 
experimental situation corresponding to a hybrid realiza- 
tion of the field configuration. We shall suppose that 
B = Bb is a static field, while E = is a light field, 
linearly polarized along the direction e, the light beam 
direction k being taken parallel to b. With this hybrid 
"B-field light-field" Hamiltonian, one can satisfy exactly 
all the conditions required in ref. [l|, by allowing e to 
rotate around 6 || A: at the angular speed a. This hybrid 
Hamiltonian 

H,,yhiF,t)^-fFnBF-b+—gUF-e)^, (9) 

is identical to the Hamiltonian of Eq.(l), if one makes 
the correspondences 



and 



A ^ ^gl/ipB. (10) 



We set 7fB = gp^^B, where jsB is the Larmor angular 
frequency of the electron and gp = 2(F — /)/(2/-|-l). (Of 
course the fact that the internal angular momentum F, 
relative to a given hyperfine sub-level, can be treated as 
an isolated spin implies the reasonable assumption that, 
for values of A ~ 1 , the Larmor frequency is much smaller 
than hyperfine splitting.) The use of an ac-light field in- 
stead of a static electric field has other advantages be- 
sides the magnitude of the Stark coupling. By adjusting 
the laser detuning, it makes it possible to apply a light 
field upon a single ground state hf level, the second one 
remaining spectator and providing the phase reference. 
The expression for A also shows there are two indepen- 
dent ways of reversing the sign of the Stark versus the 
Zeeman-coupling: one can reverse either the sign of B or 
else that of the laser detuning. 



B. Solving physical problems raised by the 
instability of the "dressed" atomic ground state. 

The main drawback of the light-shift method for ob- 
taining an effective E field, is that it generates an insta- 
bility of the "dressed" ground state hf sublevel, which 
is going to simulate our isolated spin system. This in- 
stability is best understood within a fully quantized de- 
scription of the atom interacting with the radiation field. 
The vector state of the "dressed" atom of interest can 
be written as: \FnSi/2) \N euj) where N is the num- 
ber of photons with energy huj and linear polarization e 
in the coherent light beam. The light shift is associated 
with the two successive virtual electromagnetic transi- 
tions: \FnSi/2) <»\Niuj) ^ \J^nPi/2) <E)\N - liu}) =^ 



\FnS1j2) ® \N eu)). More precisely for an ap- 

propriate detuning, is returning by stimulated emission 
to the same ground state hf sublevel |F 715*1/2) but with 
its energy "light-shifted". An alternative route for the 
excited state \TnPi/2) is to make transitions to ground 
state hf sublevels, by an energy conserving spontaneous 
light emission, leading to the infinite set of final states: 
\F' nSi/2) (»\N - liuj;! i' uj') where w' 7^ w and e' ^ e. 
We see easily that all these final states are orthogonal 
to the initial "dressed" ground state hf sublevel even if 
F — F' . In conclusion, our candidate for an isolated spin 
system is unstable with a decay rate Tdec which scales as 
(ri^/A^)r„Pj^2 J where T np^/^ denotes the spontaneous 
emission rate of the excited state 7iPi/2- 

At first sight, since the light-shift involves the inverse 
of the detuning, while the decay rate scales as the square 
of this quantity, the choice of large detunings would seem 
the most appropriate. In fact, this is not so. Indeed, for 
?iA much larger than the hf frequency splitting of the ex- 
cited Pi/2 state, AWp, the two hf states contribute with 
nearly equal magnitudes but opposite signs, in such a 
case the light-induced quadratic spin-coupling vanishes. 
Therefore, the quadratic coupling is easier to achieve at 
detunings comparable to AWp for heavy alkali atoms 
(Rb, Cs) or odd isotopes of alkali-like ions (Ba+, Hg+), 
which have the largest hf splittings. Since the P3/2 hf 
splitting is smaller than that of P1/2 by a factor 5 in Rb, 
the light beam should be tuned preferably close to the 
Di resonance line. Detunings of the order of this split- 
ting (0.816 GHz for ^''Rb) appear to lead to a reason- 
able compromise solution between decay rate and light 
shift magnitude (comparable to the Zeeman coupling for 
a magnetic field in the 0.1-1 mG range, see Sec. IV). For 
a correct evaluation of A, we have to include the contri- 
butions from both hf nPi/2 levels, with their respective 
detunings Ap.jF = uj — ujp.j^. The explicit formula giv- 
ing His{F) in the case of an alkali ground state hf level 
F {gp = 1/{I + 1/2)), lightened by a linearly polarized 
light beam becomes |4ll, l42i | 



HisiF) 



A24 -Hirp/2 



i2,2 



iVp/2 



{l~gl{Y-ef) 
gU^-ef. (11) 



A remarkable feature of this expression is that, as an- 
nounced, the two hf nPi/2 levels quadrupole contribu- 
tions cancel out if the frequency detuning is much larger 
than the nPi hf splitting, A2,2 — A2,i. In the sum with 
equal weigths of the hf levels numerators the nuclear spin 
dependence disappears; only remains the effective tran- 
sition dipole cx s. The instability discussed above has 
been accounted for by adding to the energy denominator 
+iTp/2, see ref.fH]. 

Ignoring for a moment the natural width compared to 
the detuning, the quadrupole coupling term is given by 
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an expression identical to Eq. (8). With the definitions: 

Wis=fi^gl and 1/A = I/A2.2 - 1/A2,i. (12) 

we recover the expression of the "B-field light-shift" 
Hamiltonian; 



hyb 



h^F BF - b+WlsiF ■ 



(13) 



deduced from Eq.® by substituting A in place of A. 

Let us now turn to the physical effects induced by 
the imaginary part of the energy denominators. The 
most dangerous effect involves the imaginary part of the 
quadrupolar coupling given by the expression 



^2,2 



up to third order in Fp/A. This contribution modifies 
during the quantum cycle the structure of the atomic 
wave function, which is a linear combination of m- 
dependent angular momentum states, and hence this in- 
validates the Berry's phase derivation of ref. [l[. There 
is, fortunately, a remedy to this problem. It is to tune 
the laser frequency midway between resonance with the 
two hf states as illustrated in the insert of Fig. ([31). Then, 
the imaginary part of A is vanishing. From the remaining 
scalar part results a decay rate F^ec of the dressed atoms, 

A^ 1 = A^ 2 implying |A2^2| = |A2^i| = 27rAWp/2,(14) 



and Wi^ — 



27rAyyp ' 



^dec — 



(27rAyyp) 



(15) 



From now on, therefore, we shall assume that the laser 
beam frequency detuning satisfies condition (I14[) and we 
rely on Eqs. (|T3|) and (|15l) . to suggest experiments on al- 
kalis. We arc left with an "isolated" spin with decay rate 
Tdec- For dressed ^'''Rb atoms, the spectroscopic param- 
eter governing the instability turns out to be small: 



laser beam, with ellipticity axes at 45° of the x and y 
axes, i.e. u = 7r/4, and a B-field along the beam: 



H{Wis,B,6e)/n = Wis [ ^sm{26e){F,Fy} - ^F^ 



-{Wis cos(2,5e) + jfB)F, + Wis ^,(17) 



with Wis 



2ttAWp 



(18) 



C. The special case of ^^Cr chromium isotopes 
with S=3 

Chromium atoms in their ground state have also L = 
like alkalis, but they possess a large electronic spin 5 = 3 
and an isotope of 84% natural abundace, ^^Cr, without 
nuclear spin. The ground state ^5*3 is coupled by dipole 
transitions to different P-states (^^2,3,4) which have been 
used for cooling the atoms in an optical dipole trap [sl] . 
With linearly polarized blue lasers suitably detuned from 
these transitions it is possible to induce Stark shifts pro- 
portional to in the ground state. For instance in the 
particular case of a small detuning with respect to the 
J conserving transition ^6*3 ^ ''P3, equation ^ could 
be adapted by replacing F by J and performing some 
angular momentum algebra. Like in the case of alkalis, 
detuning and light power have to be adjusted to mini- 
mize the instability of the dressed atomic state. How- 
ever, as in chromium the fine structure of the P-states 
being about 300 times larger than the hf structure of the 
alkali P-states, those restrictions are much easier to sat- 
isfy. Importantly, it is also much easier to prevent the 
quadratic coupling effect from being affected by the in- 
stability. This quadratic Stark shift has been observed, 
and exploited for the study of spin-3 Bose-Einstein con- 
densates in the geometry E 1 1 B [3^] (where it was termed 
"quadratic Zeeman effect" ) , but we are interested here in 
the geometric phases in the configuration E _L B. 



27rAyVp 



= 7 X 10 



-3 



(16) 



instead of 6 and 5 xlO"^ for ^Li and ^'^Na. As shown in 
Sec. IV, *^Rb thus provides a unique possibility to test 
Berry's cycles for 5 = 2 in the range < A < 1. 

In the last section we shall be interested in the case of 
a spin-1 performing a quantum cycle not in the Hamil- 
tonian parameter space but in the density matrix space 
described by the atomic polarization and alignment. The 
hybrid light-field B-field Hamiltonian considered in this 
section is adapted to realize a cycle of this kind. The 
laser beam in this case will be equally detuned from the 
two hf lines starting from the F — 1 ground state, there- 
fore A^2 = ^ii. It is easily verified that all preceding 
equations remain valid if A12 is changed into A21 and 
A22 into All. We shall need the explicit hybrid Hamil- 
tonian for _F = 1, in the case of an elliptically polarized 



IV. POSSIBLE ATOMIC INTERFEROMETER 
MEASUREMENT OF BERRY'S PHASES FOR 
SPINS WITH QUADRUPOLE COUPLING 

Among several possible tests of the most striking re- 
sults associated with the quadratic Stark coupling, there 
is an important goal which is to first observe the Berry's 
phase generated by the rotation of E around B. 



/3(m = 0,A) = - (pp{0,X)da 



(19) 



This is one of the remarkable effects signalled out in Sec. 
II, that we expect to be large for 5' = 2 when A is close 
to 1, as highlighted by figure 1. The effect is particularly 
noteworthy for a spin larger than one, initially in an \S, 0) 
substate, i.e. having its quantum averaged polarization 



11 



along B cancelling out before an adiabatic cycle starts 
and after it stops. Therefore, the purpose of this section 
is to make precise suggestions for a measurement involv- 
ing the rotation of E around the orthogonal B field, when 
both couplings, linear and quadratic, are of comparable 
magnitudes. 



A. Experimental compromises for spin-2 
measurements with matter-wave interferometers. 

Measurement of the phase difference acquired during 
the evolution of a quantum state requires a phase ref- 
erence. Therefore the methods of matter- waves inter- 
ferometry, first employed with neutron beams (i^ . |45| . 
then adapted to magnetic resonance [111, and now 
the subiect of outstanding developments in cold atom 
physics 131, are especially well suited to observe Berry's 
phases. Atomic interferometry has been used before to 
measure a topological phase p8l - l52j , but not in the case 
of a quadratic spin coupling, where the parameter space 
cannot be reduced to the surface of a sphere. 

The initial atomic state is represented by a coherent su- 
perposition of two states having different energies, such 
as for instance two hf substates of an alkali atom. The 

main difficulty arises from the decay rate Tdec = itj— Tp 
"2,1 

(see section III.B Eq. p^ V For concreteness, we shall il- 
lustrate our proposal on the case of the two hf substates 
of ^''Rb, |F = 1, TO = 0) and \F = 2, to = 0). These have 
been chosen because first, only one of these two mixed 
substates acquires a Berry's phase; secondly, as we shall 
show, the range of parameters leading to A w 1 looks 
achievable for this isotope. By making a specific choice 
of the experimental parameters, we have found that an 
interferometer measurement of the peculiar features of 
the 5 = 2, TO = Berry's phase is feasible within the 
present state of the art. However, as we shall see, there 
is actually little freedom for organizing the quantum cy- 
cle, if one wants to perform a measurement, free of non- 
adiabatic corrections. 

i) Instability versus adiabatic requirements 
The relevant range of the A parameters (0 < A < 1-5), 
has to be reached with moderate laser intensities, in order 
to limit the decay rate of the dressed atoms. It requires 
the magnetic field to be small (~ 1 mG), but not much 
smaller, in order to keep the field homogeneous over the 
atomic sample. The problem of the dressed atom losses 
would favour rapid quantum cycles but they should not 
spoil the validity the adiabatic approximation, in par- 
ticular during the preparation of the coherent quantum 
state before it is submitted to the rotating electric field. 
It is then crucial to tame out the non-adiabatic oscilla- 
tions generated by a linear ramping of X{t) from zero to 
a value ~ 1. An efhcient remedy, given in is to use 
a Blackman pulse-shape for the time derivative X{t) (See 
II.D) 

As regards to the non-adiabatic correction associated 
with the angular velocity of the E field rotation, they 



are governed by the convenient parameter 77 = a/j^B. 
We have shown in our previous work [1] that all the cor- 
rections to the Berry's phase of the order 7y2"+i vanish 
for n > 1, if X{t) is given by the magic value X*{2,ri{t)) 
( See Section IID.) We shall see that the set of values 
77 = 0.3 and A* (77) — 0.830, is an element of the overall 
compromise to be elaborated in the next subsection, 
ii) Absence of quadrupolar light shift in the F = l,m = 

quantum state. 

The cancellation of the quadratic light shift in the |1, 0) 
quantum state of the ^^Rb atom results from its absence 
of TO = ±2 sublevels. This simplifies the preparation of 
the coherent state which becomes a superposition of this 
unperturbed substate with the state perturbed by the 
laser field. This latter can be written |4'(2,0;t = 0)) = 
|7/'(2, 0; A(0))) in the initial state of the quantum cycle; 
and is thus an eigenstate of H{'B{t),'E{t)) defined by Eq. 
([T^. (Note that it is now useful to specify the spin value, 
since we are dealing with an admixture of two different hf 
states, hence belonging to different internal spin spaces), 
iii) Parameter compromise for a simplified S = 2, m = 
Berry's cycle performed on ^'^ Rh: one example. 

Up to the end of this subsection and only here, we shall 
analyse the simplified S = 2, to = Berry's cycle defined 
by the boundary conditions: A(0) = A(Tc) and ci{Tc) — 
q(0) — TT and made the further assumptions A(t) = and 
ait) = tt/Tc. 

Table I presents the relevant parameters for ^^Rb. For 
A2,i = iAWp/?i = 27r X 0.408 GHz and = 27r x 
0.68 MHz corresponding to a laser intensity of « 0.37 
mW/cm^, namely 0.113 times the "saturation" intensity, 

1 sat defined such that_ri^^( = rp/8; the quadrupolar 
Stark coupling, il^gln/A, amounts to 27rx 575 Hz, while 
the linear Zeeman coupling is 27rx 700 Hz/mG. In order 
to explore the interesting domain 0.1 < A ^ 1.5, the 
applied magnetic field should lie in the range 10 > B > 
0.65 mG. 

For constant angular velocity, the time needed to per- 
form one quantum cycle is = 7r/d = -k /rj^pB. If we 
want to keep the signal loss per cycle arising from the 
light- induced Stark coupling smaller than 1, the condi- 
tion to be fulfiled is TdecTc ^ 1, or 



4n- ( 



77 '2ttAWp 



) < 1, (20) 



using Eq. ljlOl and I15 p . Although the laser intensity is 
involved in both the expressions for and T^ec, H disap- 
pears from TdecTc^ but the intensity selected determines 
the magnetic field range and the minimum time needed 
for a measurement, which are interrelated together: the 
higher the field, the shorter the time. As noted before 
^^Rb is among alkali atoms the most favourable one with 
Vp /2tt AyVp = 7 X 10~^. Therefore, one can only select 
the value of the ratio A/t;, preferably not exceeding a few 
units, for condition (pil)) to be satisfied. If one wishes to 
take advantage of the "magic" value property of A, one 
is led to choosing a rotation speed of the E field moder- 
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TABLE I: Experimental parameters for a Berry's phase measurement with *^Rb atoms supposing A2,i — — A2,2, and using 
time- independent "magic" conditions -q = 0.3, A* (77) = 0.83030. Values are given for two different laser intensities / expressed 
in terms of the "saturation" intensity Isat, defined as O^at ~ Tp/S. For the Di line of Rb Isat ~ 3.3 mW/cm^. The product 
TdecTc does not depend on / (see the text and Ea. (|20p ). 
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FIG. 3: (Color online) Interferometric measurement of the Berry's phase for a spin 2 (m = 0), chronology of the interferometric 
cycle including the quantum cycles, time in arbitrary units. The two phase-coherent Ramsey pulses are resonant for the 
1,0 — >■ 2,0 transition of unperturbed *^Rb atoms. The quantum cycle of atoms in the state |2,0) starts with a ramping-up 
of the intensity of the off-resonant laser field, then, the rotation of its linear polarization starts smoothly towards a maximum 
angular speed. The whole operation is followed by the time reversed one to return to the initial state. The temporal dependences 
of both the light shift - normalized by the Zeeman shift, X{t)- and the angular velocity of the polarization rotation - normalized 
by the Larmor precession angular frequency, ri{t) - as well as their relative magnitudes are the parameters which govern the 
non-adiabatic corrections to the Berry's phase in a critical way. By using Blackman pulses for both X{t) and rj{t) ("oscillation- 
taming") and by satisfying the magic relation \{t) = X*{ri{t)), we show (see the text) how it is possible to keep these corrections 
well below the 0.1% level, with 40% of the initial cold-atom cloud surviving the light-induced decay. Insert: relevant atomic 
levels. 



ately large, (77 = 0.3, X*{r])/v = 2.76 for ^^'^Rb). Finally, 
this choice of experimental conditions (see Table I) leads 
to an acceptable signal loss for a single quantum cycle, 
~ 24%. 



B. Towards an empirical determination of Berry's 
phases, free of non-adiabatic corrections to the few 
0.1 % level 



1. Organization of the interferometric cycle 



In the next subsection, we propose a realistic precise 
timing for the quantum cycle including now a discussion 
of the atomic loss during ramping up and down of the X{t) 
parameter. This operation is of a critical importance to 
avoid non-adiabatic corrections. 



The interferometric cycle can be organized according 
to the method developed for high accuracy measurements 
with hfs measurements in cold atom atomic clocks to pre- 
pare the unperturbed coherent state oc (|1, 0) -I- c |2, 0)) 
{e.g. [5^). Starting from the pure |1,0) state in pres- 
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'W2,2(t) 




FIG. 4: (Color online) Time dependence of the diagonal {S — 
2, m = 2) matrix element of the Hamiltonian in the rotating 
frame, during the whole quantum cycle (time unit: (7s_B)~^). 
Black curve: linear ramping for A(i) and square pulse for 
a{t). Red curve: taming of the non-adiabatic oscillations 
(OT procedure) by using Blackman pulses for both \{t) and 
a{t). 

ence of B and without the laser beam, hence vi^ithout E, 
one can use the Ramsey method [5] to prepare the un- 
perturbed coherent state cx (|1,0) + c |2,0)), i.e. apply 
suddenly a short pulse (using either a rf field resonant for 
the 1 1,0) — 1 2,0) frequency or two Raman pulses [s^j). 
After a certain delay (time of free evolution, T, long com- 
pared to the durations of the rf pulses and the quantum 
cycle), the quantum cycle starts. The laser field respon- 
sible for E is applied progressively and follovifed later on 
by the polarization rotation. It lasts for a duration Tc. 
Then, detection takes place. This consists in the applica- 
tion of a second Ramsey pulse, in presence of the B field 
and without E. The Berry's phase is obtained from the 
probability of finding the atoms in the |1,0) (and |2,0)) 
state at the end of this pulse. The latter is measured 
by fluorescence detection following resonant excitation of 

the 5S'i/2,F=i - 5P3/2,F=o (and 5S'i/2,f=2 - 5P3/2,f=3) 
transitions. The cycle is depicted in Fig. [3] 



2. Implementing the quantum cycle via 
an oscillation-taming procedure 

For clarity we suggest splitting the quantum cycle into 
three steps in which only one parameter at a time is 
varied. During the first and the last steps, A alone is 
varied from to Aq, while during the second step a is 
varied from to tt. Moreover, we decide to comply with 
the recommendation given and justified in [l[ to satisfy 
the adiabatic approximation, for preparing the quantum 
state ^'(2, 0; Aq) at the beginning of step ii) with its po- 
larization (Sz) very close to its desired adiabatic value, 
p(2 0;Ao). We avoid the discontinuities in the variation 
of X{t) by suitable tailoring of its shape: A(t) will be as- 
sumed to be described by the primitive of a Blackman 
pulse. During step 2) the time variation of a is supposed 
to be represented by the Blackman function /(s) (Eq. [5]); 
this is what we have termed the oscillation taming (OT)- 
procedure. In addition, it is advantageous and possible 



to adjust \{t) so that at any time it coincides with its 
"magic value" associated with the angular velocity 77 (i) 
at that time. This means that during steps i) and iii) 
Ao is chosen equal to A'^(O), while during step ii), instead 
of keeping A constant, we make a fine tuning of A in or- 
der to satisfy Eq.(|31). In practice, the variation of A*(r/) 
in the range < ry < 0.3 is only half a percent of Aq. 
The advantage of this strategy is the complete suppres- 
sion of non-adiabatic corrections to the Berry's phase, 
which are odd under reversal of a (Sec.V.A), the even 
ones being suppressed by subtracting measurements for 
mirror-image cycles {a(t) — > —a(t)). 

Table II summarizes the chronology of Berry's cycle 
organized in three steps. For each step we give the ex- 
pression the Hamiltonian ^.{t) acting within the frame 
rotating around the z axis with the angular velocity d, 
together with the time dependence of the parameters X{t) 
and ri{t). 



3. Quantitative predictions 

Choosing Tc = 3T which implies for the angular speed 
at its maximum j^™"^ — 0.27 ( quite close to the value 
considered in Table I), we have simulated the exact time 
evolution. We have solved numerically the Schrodinger 
equation for the Hamiltonian represented by the sum of 
the three time dependent operators defined on the three 
successive time intervals. (Rounding-off theta functions 
are used to avoid discontinuities of higher order deriva- 
tives at the passage between two time-intervals.) The ini- 
tial state at i = is the pure state |2, 0). If our attempt 
to create the conditions for the adiabatic approximation 
is working, we expect {Sz) to be very close to its adiabatic 
value at the end of step i) and at the end of step ii) . The 
calculation yields a difference of only 1 part in 10"^. This 
indicates that, as anticipated from the discussion given 
m [l| (Sec.V.C), our implementation of the Berry's cycle 
leads to very small deviations from the strict adiabatic 
evolution, concerning the quantum states. 

We have also probed the ability of this implementation 
to reproduce the adiabatic approximation as regards to 
the phase. We have performed two kinds of test. First 
we have calculated the phases of the final state for two 
mirror-image cycles (d(t) — )■ — Q;(t)), still by solving nu- 
merically the Schrodinger equation, and compared their 
difference A$f) to the adiabatic Berry's phase, extracted 
by numerical evaluation of expressions ([3]). Our result, 

sin (A$z5 - (/?+(2, 0; A*) - ;3~(2, 0; A*))/2) = -0.000035, 

(21) 

is a good confirmation that, with the chosen timing, the 
corrections to the adiabatic approximation can be made 
exceedingly small. However, such an identity provides a 
determination of the Berry's phase only modulo tt. Using 
the symmetry l3+{2, 0; A*) = -,3" (2, 0; A*), the result of 
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TABLE II: Chronology of Berry's cycle organized in three equal steps using 1/jsB as time unit. The function h{t) is the prim- 
itive of the Blackman function /(s) (Eq.©), varying between and 1 over the unit time interval, h{t) — f{s)ds/ J^^ f{s)ds. 

Time interval H{t) \{t) ■q{t) = a{t) /-y^B 

step i) < t < r S, + \{t)Sl A*(0) h{t/T) 

stepii) T<t<2T {l-ri{t))S. + \*{r,{t))S^, A*(77(t)) Eq.© _^ /, (1^) 

step iii) 2T<t<Tc^3T S, + \{t)S^ A*(0) h {{n - t)/T) 



our calculation actually is: 

iA$D =/3 + 7r-3.5x 10"^ (22) 

The presence of tt should not be considered as a surprise 
since our determination of the phase at the end of each 
cycle is obtained from the argument of the wave function, 
and is therefore defined modulo In. 

To remove the resulting ambiguity of tt on the half- 
difference between mirror-image cycles, we have per- 
formed a second test. We compute the difference of the 
adiabatic phases, A^acUab directly from the exact expres- 
sion of the instantaneous eigenenergies of the |2 0; A(t)) 
state during the selected cycle, by doing well-defined 
quadratures. After comparison of both evaluations, the 
result 

^A^D - ^A$,d,,6 + TT - 2. X 10-5 (23) 

merits two important remarks. First, the implementation 
of the quantum cycle that we have selected, succeeds to 
give an excellent control not only of the quantum state 
but also of the adiabatic phase. Second, the tt incre- 
ment appearing in our numerical evaluation of Eq. (j22p 
confirms the fact that our theoretical evaluation of Acj)!) 
is defined only modulo 27r, as will also be the case for the 
experimental determination, A$^^. But, at the same 
time, we obtain the means to remove the modulo tt am- 
biguity in the determination of /3: it is enough to look 
at the difference i(A<i>^^ — A^adiab) evaluated with a 
precision comparable to the experimental one. If, as we 
expect, this is found equal to an integer times tt with 
enough accuracy (depending on the experimental preci- 
sion, but better than 10 %), then the modulo tt ambiguity 
is suppressed from iA<i>^^. 

Finally, the resulting empirical determination of the 
Berry's phase can be made free from systematic uncer- 
tainties caused by deviations from the adiabatic approx- 
imation, within an accuracy even beyond the 0.1% level. 
If we suppose that the time variation of A is realized by 
adjusting the laser intensity while keeping the B field 
constant ~ 1 mG, the whole cycle duration amounts to 
~ 15 Larmor periods, i.e. 21 ms. Obviously, there is a 



compromise between accuracy and duration: a faster A- 
ramping could be employed but at the expense of lower 
accuracy. 

4. Interferometric detection of the Berry's phase 

We can now present an explicit evaluation of the opti- 
cal signal detected for a realistic choice of the parameters 
involved in the interferometric measurement. The Hamil- 
tonian associated with the first rf pulse is written in the 
laboratory frame : 

^rai-^ + ^^i(-e™* + r+e— )• (24) 

In this subsection we shall use a unit system such that 
h — 1 and we set w = 2ttAW. The symbols ti,T2,T3 
stand for the familiar Pauli matrices and t± = (ri ± 
1x2) /2. With these notations lui, represents the cou- 
pling of the rf frequency field Brfzsinwt with the mag- 
netic dipole transition operator from the lower state 
= l,m = 0) to the upper state \F = 2,m = 0), 
and the mixed states are described by a two-component 
spinor ^'(t) = {X2{t), Xi{t)}. A straightforward compu- 
tation performed within the frame rotating at the angular 
velocity w leads at the end of the first rf pulse {t — tp^ ) 
to the following state vector : 

^,„,(^,J = {^sin(^^),cos(^^^)}. (25) 

The value of uJitpi will be specified later on. During the 
free evolution of duration T, the rf field is decoupled the 
rotating frame Hamiltonian vanishes, so that ^rot{ipi + 
T) = ^'rot(ipi), when the quantum cycle starts. During 
the quantum cycle the spinor component X2{t) acquires 
the phase <I>Jq 

'i>2%-*D + /3*(2 0;A*), (26) 

the ± index refering to the sign of 77. Simultaneously, it 
is also affected by the light- induced decay, but Xi is not, 
so we obtain 

{, sin(^:^) exp (-frfec T,/2 + ^$20) , cos(^^)}, (27) 
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where Tdec represents the decay rate averaged over the 
whole quantum cycle. 

Up to now, ujitp-^ has been considered as a free pa- 
rameter. It turns out to be advantageous to adjust it to 
make the two spinor components of equal magnitudes at 
the end of the quantum cycle. For this purpose, the con- 
dition to be satisfied is tan {uJitp^/2) — exp (rdecTc/2^ — 
1.953 in the present example, leading to arctana;itpj/2 = 
1.0975 rad, and a common magnitude of 0.455, instead 
of 0.707 when there is no decay. The signal loss resulting 
from the dressed atom instability is thus a factor 2.5. 

A second 7r/2 Ramsey pulse, can be applied to the 
state vector and be immediately followed by the detection 
process. The measured quantities are the probabilities Vi 
and 7^2 of finding the atom in states F = 1 and F — 2 
respectively 



(1-COS$20'^) 
(1 -^COS$20'^)- 



(28) 
(29) 



After making the two measurements for 77 > and 77 < 
and using Eq. I^E^, one can extract 2/3(2 0; A*) modulo 
2'K. We recall that the ambiguity of tt appearing in this 
determination of /3(2 0; A*) can be removed by combining 
this result with a calculation based on the knowledge of 
the eigenenergies. We see that, at the price of reducing 
the number of detected atoms by a factor 2.5, the fringe 
visibility can be kept very close to unity. 

Subsections A and B mainly refered to cold alkali 
atoms or trapped alkali ions, offering integer spin val- 
ues between 1 and 4. However, in the expanding family 
of laser cooled and trapped atoms there is the much less 
familiar case of chromium which, we believe deserves a 
special attention. ( See section III.C). 



V. BERRY'S PHASES FOR SPINS 3/2. 
POSSIBLE MEASUREMENTS 

A. Theoretical background 

Although half-integer spins have no fully symmetric 
state like the m = substate of integer spins, their 
Berry's phases for quadratic Hamiltonians are interest- 
ing in their own right. With S=3/2 we have the lowest 
spin state for which the dimension of the density ma- 
trix space exceeds that of the parameter space. In or- 
der to obtain a spin Hamiltonian able to generate the 
Aharonov-Anandan phase, one would have to include an 
octupole spin coupling. On the other hand, analytical 
Berry's phase formulas relative to the S = 3 /2 q uadratic 
spin Hamiltonian have already been derived [29| . In Ap- 
pendix B, we give explicit expressions for f (m. A) and the 
polarization p{m, A) obtained by rewriting the results of 
ref [29] within the notations of [l| and the present paper. 
The results are displayed in Fig. 5. 



The most remarkable feature apparent in this figure 
concerns the level m = — i. The polarization p(— 1/2, A), 
governing the size of Berry's phase, becomes > ^ for A > 
2 indicating a strong mixing with the state m = | . At the 
same time, this level remains at a distance ~ 1 {jsB unit) 
from all other levels. This situation looks favourable for 
an empirical determination of Berry's phase /3(|,— ^). 
Moreover, /3(|, — i) differs widely from the spin one half 
Berry's phase /3(^, — ^). This difference provides the pos- 
sibility of achieving a maximum entanglement of three 
non- correlated spins by the method described in section 
VI. of reference [l[ . More details about this "holonomic" 
three Qbits entanglement are given in Appendix A. 



B. Possible measurements 

There are two kinds of 3/2-spin systems for which 
measurements of Berry's phase might be considered : i) 
the '^^Cl nuclei embedded inside a monocrystal matrix 
21 1, [131 and ii) the ^°^Hg mercury isotope in a vapor 
26| or a cold atom optical trap 53] (which could apply 
as well to odd isotopes of alkali-earths atoms like ^^^^Ba 
and ^■^^Ba placed in an optical dipole trap). Up to now in 
those systems, none of the experiments performed so far 
corresponds to what should be ideally realized for testing 
original features of Berry's phase. 

The first experimental approach corresponds to the 
nuclear quadrupole resonance (NQR) of ^^Cl nuclei (1= 
3/2, fi = 0.82 /i„) of an oriented axially-symmetric sin- 
gle crystal of sodium chlorate put inside a sample ro- 
tor [2i|. In this case, the quadratic spin coupling re- 
sults from the interaction of the '^^Cl nuclei quadrupole 
{Q = — O.OSx 10~^^cm^) with the local axially-symmetric 
electric field gradient. However, up to now, measure- 
ments have been performed only without the linear cou- 
pling needed to lift the level-degeneracy. To investigate 
the interesting range A ~ 1 the magnetic field B (about 
10 G) should be applied along the direction z of the 
first rotor axis (ssj , the quadrupolar axis of the rotating 
monocristal being oriented perpendicularly to B. The 
second rotor permits the precession of B about the ref- 
erence axis z. This set-up would be adequate for precise 
verifications of the Berry's phase and its non-adiabatic 
corrections. 

In the second situation relative to ^"^Hg, the conditions 
corresponding to the hybrid "B-field laser-field" situation 
satisfying the E • B = hypothesis, have been achieved 
in optical pumping experiments of early 1970's [2^. The 
quadratic spin coupling was adjusted as wanted over a 
large range of coupling ratios A, (although there was no 
laser available at that time) but, though the applied fields 
were time-dependent, they were far from being applied 
adiabatically. 

Measuring the Berry's phase could be performed by 
preparing ^'^^Hg atoms in a coherent superposition of 
two Zeeman substates, the most interesting one being 
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S=3/2 

£'(m,X) p(m,A) 




FIG. 5: (Color online) Reduced energies £{m,X) and polarization, p(m,\) for S — 3/2 within the interval < A < 4.( The 
results for A < are obtained by reflexion about the axes origin). The findings relative to the level m = —1/2 for 2 < A < 4 
are quite remarkable. This level remains far away from the other ones and at the same times its polarization p(— 1/2, A), giving 
the size of the Berry's phase, reaches values > 1. This indicates a strong mixing with the level m = 3/2. 



'^coh = ^|3/2,-3/2) + |3/2,-l/2) for A > 0. To this 
end one can start from the pure |F, mp) — |3/2,— 3/2) 
substate prepared by optical pumping and apply either 
a rf-frequency field or a pulsed modulated light beam to 
perform transverse optical pumping [26]. The fictitious 
electric field can be created using adiabatic application 
of a laser beam detuned from the hyperfine lines of the 
"'^iS'o — ^ ^Pi transition. The resulting ground-state insta- 
bility will be mild compared with the case of ®^Rb for 
two reasons. First, the electron-Zeeman coupling is re- 
placed by the much smaller nuclear-Zeeman coupling, so 
that for a given detuning, much weaker radiation fields 
are needed to realize A « 1. Secondly, the life time of 
the ^Pi Hg excited state is longer and the hfs splitting 
larger, allowing one to select a frequency detuning of a 
few gigaherz, the ratio Tp/AWp is thus reduced by one 
hundred. Finally with all these parameters combining 
in a favourable way, in Hg the signal loss per cycle (Eq. 
pO|) becomes unsignificant. 

After one adiabatic quantum cycle, the phase shift 
induced by the time-dependence of E, expected to be 
n„7r(l -p(-3/2,A) -I- p(-l/2, A)) (Eq.®), is close to 
TT over the range 2 < A < 4. There are several rel- 
evant signals signals: the modulation of a transmitted 
resonant probe pulse of low intensity at the frequency 
7_FB/i-i(f (-1/2, A) _ £:(_3/2, A)) characteristic of the 
dressed atomic coherent state or the optical rotation of 
a linearly polarized beam tuned off-resonance of one hf 
component of the ^Sq—^ Pi transition. If the probe pulse 
is applied at the end of the quantum cycle, the quan- 
tity to be measured is the phase shift of this modulation 



generated by the rotation of the E-field. Since the sign 
of this shift changes with the direction of the rotation, 
making two consecutive measurements, with opposite ro- 
tation velocities, should yield the Berry's phase. When 
both results are combined, the modulo tt ambiguity has 
to be resolved in the same way as indicated for the *^Rb 
interferometry experiment (cf. Sec. IV). 

A different measurement scheme could exploit optical 
pumping with polarization modulated light, which can 
generate selectively high-order coherences [5a, [53] ■ 



VI. ANANDAN-AHARONOV GEOMETRIC 
PHASE FOR S=l 

In this section we shall deal with the AA geomet- 
rical phase, still unobserved in the case of a spin-1. 
This phase is generated by quantum cycles along closed 
circuits drawn upon the space E(p) of the pure state 
density matrices p{t) = |^(t))(5'(i)|. During the cy- 
cle, at any instant t, the parallel transport condition: 
(\I'(i)|^^'(t)) = has to be exactly satisfied. After con- 
sidering the general form of the parallel transport Hamil- 
tonian 28:], we shall show that the hybrid B-field light- 
shift Hamiltonian can be tailored to suit to this form by 
adjusting the time-dependences of its parameters. This 
allows us to suggest a method for performing the mea- 
surement of the A A phase for spin-1, valid for example 
in the case of ^^Rb atoms in the F = 1 hf state. 
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A. Aharonov-Anandan's versus Berry's phase as 
physical objects 

In the spin-1 case the density matrix is com- 
pletely determined by the knowledge of the polar- 
ization vector Tr{p{t)S)/h and the "alignement" ten- 
sor Aij = Tr({Si, Sj})/h'^ . By performing an ap- 
propriate rotation upon the spin system R(t) — 
R{z,(f{t)).R{y,9{t)).R{z,a{t)) involving the three Euler 
angles, the tensor Aij can be put under a diagonal form: 
A= z z + ^{1 + sin X X + ^{1 — sin C)y(>^y with the 
polarization lying along the z axis, p = cos(C) z where 
— 7r/2 < C < '^/^- The angle together with the Euler 
angles, provide a system of coordinates for which, 
by construction, is isomorphic to CP^ . In Reference [1^ , 
G. Gibbons and one of us (CB), we have derived the ex- 
plicit expression of the AA geometrical phase in terms of 
the above set of coordinates: 

= ^ cos(C) cos(6i) - (1 - cos(C)) da. (30) 

The closed loop C is specified by the following constraints 
within the time interval < t < T^: ({Tc) = C(0), -f < 
at) < f , e{n) = 0(0), < 9{t) < n and ^(T,) = (^(0) + 
2n^'!T, a{Tc) = a{0) + UaTr, where and Ua are non- 
vanishing integers. By setting A = 2 tan(C) and using 
the method described in our previous work [l'|, one finds 
easily that Berry's phase for S — tti — 1: /3(1,1) = 
Paa mod{2Tr). 

We would like to analyse the respective behaviour of 
/3(1, 1) and 13 a A in the limit ( = ^ + 0^ or equivalently 
A — > +0O. It is easily seen that our set of coordinates for 
CP^ is singular for C = since p is vanishing and the 
alignement tensor A reduces to z(^z + x(>^x. Within such 
a density matrix configuration it is not possible to define 
the Euler angle a; this is like the longitude for spherical 
maps which is undefined at the north pole. However, 
the problem disappears when = ^ — e, e being a small 
positive real, so that Paa, which is designed to be free 
of non-adiabatic corrections, is well-defined by the closed 
loop integral (PU)) . provided the conditions listed above 
be satisfied. 

The case of Berry's phase is more delicate, since lev- 
els S' = m = 1 and S* = 1 , m = become degenerate 
in the limit A — > +oo and are already very close for 
A > 2. As a typical example, let us consider the non- 
adiabatic correction A/3 given by equation (73) of section 
V.C of T|. It can be easily calculated in the present case: 
A/3 = (cos6i(^ -h d)p(2)(i,A) with p'-^Hl,X) = 

A -I- C'(l/A'^). The parameter /i, governing the non- 
adiabatic eff'ect is given by: = — sin 95/(73 i?). One 
sees clearly that the non adiabatic correction A/3 is lit- 
erally exploding when A — >■ +00. Quite remarkably, the 
situation is very different in the limit A — )• — 00: the sep- 
aration between the levels S = m = I and 5 = 1 , m = 
grows like —A and p^^^l. A) = O (l/A'^), so the non- adi- 
abatic corrections can be ignored. Concerning the AA 
phase, there is practically no difference between the two 



limits C ±7r/2. The above results suggest that al- 
though /3(1,1) and Paa are identical mathematical ob- 
jects, their actual measurement will raise very different 
physical problems, as this will be confirmed in the fol- 
lowing subsections. 

B. The spin-1 "parallel transport" Hamiltonian 

In reference ^28'], we have constructed a Hamiltonian 
-ffy (t) which performs exactly a parallel transport around 
the closed circuit C drawn upon CP^. H\\{t) must 
satisfy at any time t, the parallel transport condition: 
(^'(t)|i?l|(i)*(t)) EE Tr{p{t)H\^{t)) = 0. In the present 
section, we shall limit ourselves, for the sake of simplicity, 
to closed circuits where a{t) and C,{t) are the sole time 
dependent CP^ coordinates. To proceed, it is conve- 
nient to introduce the "rotated " Hamiltonian: H\\ (t) = 
U^R{t))Hiiit)U{Rit)) with U{Rit)) = exp(-i5^a(t)). 
In reference [28[, using equations (49) and (A5) it was 
found that (t) takes the simple form: 

^11 (t) = ^{S., Sy} - ^ cos{C{t))S'^ + a{t)S.. (31) 

In the present context we are going to consider i/y (t) 
just as an "ansatz" and show that it does possess all the 
desired properties. To this end, we shall need the "rotat- 
ing" frame Hamiltonian, which governs the evolution of 

$(t) ^ uHRimity. 

Hiiit) - j^CmS.,Sy} - ^cos(C(t))5,^ (32) 

From the identity: {S^j: + iSyy - {S^-iSy)'^ = 2i{Sx,Sy}, 
it follows that i?||(t) does not mix the state |l,m = 0) 
with the states \l,m = ±1). As a consequence, if ^'(i) 
satisfies the initial condition ^'(0) — |1, 0), it can be writ- 
ten as a two-component wave function (C(l), C(— 1)). 
Using standard text book formulas, it is then easily found 
that to derive its evolution one can replace, in H\\(t), 

the matrices {Sx, Sy} and S*^ respectively by the 
Pauli matrix ay and the 2x2 unit matrix. Exploiting the 
analogy with a Ramsey pulse, one arrives at the following 
expression for 4'(t): 

^ (cos(^), sin(^)) exp -z da cos(C). (33) 

Using equations (|3T|) and ([SS)) one gets: 
{^{t)\Hii{t)^{t) = (*(t)|i7y(i)*(t)) = 0. In other 
words, we have shown, as announced, that the Hamil- 
tonian i?y(t) = U{R{t)H\^{t)W{R{t) performs a 
parallel transport along the particular closed circuit 
C considered in this section. As a final check, let us 
calculate the AA phase in the laboratory frame. Writing 
^'(T,) = exp(-|(a(re) - a(0))^,)$(Te) and using the 
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phase shift arg(^(Tc)/\&(0)), deduced from equation 
one arrives at the final expression for the AA 

phase: 

Paa= [ 'dt(cos(C)-l)«(i), (34) 
Jo 

which does agree, as expected, with the general formula 
giving ([50)1 in the particular case (p = 0. 

C. Parallel transport with a light-shift 
Hamiltonian 

By adding to the Zeeman Hamiltonian the light-shift 
Hamiltonian His given by equation ([8]), for a beam po- 
larization e = '^(^ ~^ exp(ix)y), one gets by a simple 
calculation a possible candidate for an experimental re- 
alization of iJ|| (t) : 

Hexp = B Sz + 

+Wis (^sin(x) S, + ^ cos(x) {S,,Sy} - ^S^^ . (35) 

Making the connection with the notations of section III: 
X = 7r/2 — 2Se and S — TiF, we recognize Eq. ([T7| where 
we have dropped the c-number contribution which does 
not contribute to the geometric phase. By performing 
the identification: H^xp = H\\(t), one obtains a set of 
three relations which open the road towards experimental 
realization of the spin-one AA phases: 

sin(2,Se) = :t-^77T, 
2a cos(C) 

Wis = 2dcos(C), 

75^ = d(l - 2cos(C)cos(2^e)) • (36) 

We have found that physically satisfactory solutions of 
these equations do exist for the following simple real- 
ization of the closed circuit C, within the time interval 
< t < Tc, namely 

C{t,TcXmax) = 4C„Q:r^(l - ^), a(t,Tc) = nalTY^ 

(37) 

with < Qmax < 7r/2 and Ua an arbitrary integer. 

Figure 6 represents the three physical quantities which 
determine H\\{t)^ in terms of the reduced time variable 
t/Tc, for two typical values of Cmax, 7''/4 and 7r/3. The 
first curve displays the angular variable 5e, which spec- 
ifies the elliptical polarization of the light beam. The 
second and the third curves give in terms of the angular 
velocity d, respectively the ac Stark shift and the Larmor 
frequency associated with the external B field. Since we 
have chosen the angular velocity d time-independent the 
total magnetic field remains constant during the quantum 
cycle: the variation of the external field compensates the 



effective light-field contribution (Eq. [36]) . therefore the 
total field does not vanish during the cycle. 

One should keep in mind that the ac Stark shift intro- 
duces an instability of the the atomic ground state. For 
a particular detuning configuration, the instability of an 
alkali ground-state with 7 = 3/2 can be accounted for by 
adding to H^xp an anti-Hermitian operator proportional 
to the unit operator 1: H^xp = ^f||(i) — fr^ecl with 
4 27rAVVj. ^isi Fp and AWp being the decay rate 
and the hyperfine splitting of the P1/2 state, fEg fTS)) . It 
is clear that such an instability does not affect the phase 
shift of the surviving atoms at the end of the quantum 
cycle. Using equation ([36]) and the explicit expression of 
the A A phase ([34]) . one can calculate the average value 
of Ydec{i) over the AA cycle as follows: 

Tc{Tdec)= / dtTdec{t) = „ {(3AA+naTT). 

Jq IttAWp 

(38) 

To end this section, we quote the values of the AA phases 
when (max = 7r/4, 7r/3 and Uq, = 1 as well as the values 
of Tc {Tdec} in the case of ^^Rb, 

(max = J^I3aA = -0.4969, T, {Tdec) - 0.148, 
Cmax = J ^ /3aA = -0.8567, Te {Tdec) = 0.128.(39) 

As expected, (3aa is independent of the time scale T. The 
fact that this scaling property holds also for (F^jgc) 
is less evident. There is a strong contrast with the case 
of Berry's phase measurement, where much higher atom 
losses are unavoidable if one wants to keep the non- 
adiabatic corrections below the 0.1% level. 



D. Realization of the parallel transport 
Hamiltonian on *^Rb atoms 

We suggest the realization of AA quantum cycles on 
spin-1 for the case of *^Rb atoms in the F = 1, mp = 1 
hf state, when the sole time varying parameters are C 
and a. In practice this could be achieved by a well- 
orchestrated time-variation of the physical parameters 
Wis{t), B{t), Se{t) and a{t) which characterize the B- 
field light-field configuration with B colinear to the beam. 
Measurements can be done with a laser beam detuned 
midway from both hf components starting from the F = 
1 ground state (Ai^i = — Ai^2 = inAW-p/^) and with 
the magnitudes of B and Wis adjusted for giving to both 
B- and £'-couplings comparable magnitudes (see Fig. 6). 

A precise interferometric measurement of the AA 
phase seems possible, the interferometric cycle being or- 
ganized in a way very similar to the one described on 
Fig. 3. Now, the first Ramsey pulse prepares a coherent 
superposition oi F — l,mp — 1, submitted to the AA 
cycle, with the = 2,mi;' = state which serves as a 
reference. However, the features of the quantum cycle 
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FIG. 6: (Color online) Time dependence of the physical parameters entering into the definition of the light-shift Hamiltonian 
which enables one to ensure parallel transport. The beam helicity is cos25e, its intensity, given by Wia, Eq. (|15p . determines 
the magnitude of the quadrupolar Stark coupling and ■ysB + cos25e represents the linear dipole coupling, (for *^Rb, F — I 
atoms the beam has to be equally detuned from the two hf excited states). We assume the angular velocity constant and 
C = arctan{Wij2'ysB) ^ACma^^il - |.). 



itself will strongly differ, for comparable physical condi- 
tions (i.e. average field magnitude and light intensity). 
Instead of the progresssive application of the time depen- 
dent Hamiltonian required for measuring Berry's phase 
precisely, one can apply the "parallel transport" Hamilto- 
nian suddenly, without causing any alteration of the AA 
phase. The total duration of the quantum cycle can thus 
be reduded by one order of magnitude. Thus, the atom 
loss during the quantum cycle is also greatly reduced. 

To best illustrate the flexibility available in the choice 
of physical parameters allowing one to perform AA-cycles 
we end by a summary of the relations expressing how 
the cycle duration is linked to the laser intensity av- 
eraged over one cycle on one hand, and to the effec- 
tive magnetic field on the other hand (assuming the ^ 
and a time-dependences given by expressions (|57|) . with 

Cniax — 

(Wis) =2{/3aa + 7^) = 5.28M3, (40) 
^sBeffT, = n (41) 

In addition, we underline that, remarkably, the atomic 
loss per cycle, Tc (Tdec) , remains constant, whatever the 
absolute cycle duration (see equation (1381) ). 

VII. SUMMARY AND PERSPECTIVES 

The main purpose of the present paper is to suggest 
methods for measuring the Berry's and AA quantum 
phases, in realistic experimental conditions, satisfying 
the physical requirements formulated in our theoretical 
work P, [2^ [2^. The spins are supposed to be non- 
linearly coupled to time-dependent electromagnetic fields 



(possibly effective ones) involving the simultaneous con- 
tributions of a linear and a quadrupole coupling. We have 
avoided the situations leading to degenerate eigenvalues, 
from which result non-Abelian Berry's phases. Important 
simplifications in Berry's phase calculation result from 
our assumption that the two effective fields involved, elec- 
tric and magnetic, are orthogonal; as we have shown, the 
problem becomes then endowed with several symmetry 
properties which make it tractable for arbitrary spin val- 
ues. However, Berry's quantum cycles, performed in the 
Hamiltonian parameter space, require fulfilment of the 
adiabatic condition at any time of the spin evolution, 
a condition not easily satisfied in experiments on inte- 
ger spin values > 1 where the expected original features 
(Sec. II) deserve confrontation with our predictions. On 
the other hand, the AA phase is associated directly with 
a quantum cycle in the density matrix space. The so- 
called "parallel transport" condition is satisfied directly 
without the help of the the adiabatic approximation. The 
construction of the spin-1 "parallel transport" Hamilto- 
nian - described in a simple case in section VLB - is a 
rather difficult task. Instead of a fixed number of external 
physical parameters it involves 45 parameters necessary 
to define a closed circuit in the density matrix space. 
This construction was performed initially in ref.Q for 
spin 5 = i and, later on, in ref.[2^ for 5 = 1, where the 
physical parameters involved are then the polarization 
vector and the alignment tensor. Detailed measurements 
have been performed in ref.fiBj for spin 5 = ^ but no 
such investigation exists for 5 = 1. Our goal here has 
been to find how to go beyond the limitations encoun- 
tered so far by experiments. To this end we describe 
concrete experimental situations, mainly chosen in the 
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field of atomic physics. We propose 

• a realization of Berry's cycles for spins > 1/2 having 
quadrupole and dipolar couplings to external fieldd 
with non-adiabatic correction below the 0.1% level. 

• a realization of parallel transport quantum cycles 
on the spin 5 = 1 density matrix leading to a mea- 
surement of the A A phase. 

As seen in Sec. Ill, the total angular momentum of atoms 
in their ground state are good candidates for playing 
the role of isolated spins. For both Berry's and AA cy- 
cles, a convenient experimental tool for coupling the spins 
non-linearly to external fields happens to be the "B-field 
light-shift" Hamiltonian. There are several variants: the 
quadrupolar coupling is realized thanks to the ac Stark 
shift induced by the linearly polarized light field, while a 
dipolar coupling of comparable magnitude can be ensured 
either by an an external magnetic field (Eq. p^ ) or the 
circular polarization of the light field fEg. dTTl) '). There is 
one drawback: the ac Stark shift induces an instability 
of the "dressed" ground state. The more severe problem 
lies in the fact that the quadrupole to dipole magnitude 
ratio A acquires an imaginary part which invalidates the 
derivation of Berry's phase given in ref. 1]. In the case 
of alkali atoms (say *^Rb), we have found a simple rem- 
edy to get rid of this unwanted imaginary contribution: 
it is to tune the dressing beam frequency in such a way 
that the detunings A21 and A22 relative to the two tran- 
sitions F = 1; 55 1 J-" = 1, 2 ; 5Pi satisfy the simple 
relation A21 -I- A22 = 0. Concerning the instability of our 
isolated spin candidate, this puts a lower limit upon the 
duration of Berry's quantum cycle which cannot be too 
short if one wishes to keep the non-adiabatic corrections 
below the level of 0.1%. A realistic proposal to solve this 
delicate problem is given in subsection IV. B. 

To explore the still non-revealed Berry's phase prop- 
erties expected for a spin of two, remarkable in the case 
m = 0, we suggest (Sec. IV) a variant of a Ramsey in- 
terferometry experiment made on the clock transition 
of cold *^Rb atoms. Between the two Ramsey pulses 
the free evolution of a coherent superposition of the two 
m = hyperfine substates is interrupted by implementa- 
tion of one quantum cycle in the upper state. The atoms 
interact with the off-resonant laser field, whose linear po- 
larization rotates of tt around B. The cycle is organized 
in three steps i) the laser intensity is ramped up from 
to its maximum, with its time-derivative tailored to fit 
a Blackman pulse; ii) the polarization is rotated with a 
time-dependent angular speed described by a Blackman 
pulse; fine tuning the laser intensity makes it possible 
to adjust the two-coupling ratio at its "magic" value de- 
pending on the instantaneous angular speed; iii) the laser 
intensity is ramped down to zero in the time-reversed way 
of step i) . Measurement of the phase is repeated for the 
"mirror-image" cycle (opposite rotation speed), the half- 
difference is expected to provide the adiabatic Berry's 
phase. We have simulated the experiment by perform- 
ing the numerical resolution of the Shrodinger equation 



in the rotating frame, allowing us to extract the value of 
the phase accumulated under the effect of the rotating 
E-field. We have found the result predicted by Berry's 
phase expression with a deviation of only a few 10~^. We 
conclude that, within this chronology, the empirical de- 
termination of the Berry's phase can be made free from 
any systematic uncertainty caused by deviations from the 
adiabatic approximation, within an accuracy well better 
than 0.1%. In addition, the selected timing contributes to 
minimize the cycle duration, once given the magnitude of 
non-adiabatic correction tolerated for one measurement. 
It is found that the problem raised by the state insta- 
bility inherent to the light-induced quadratic coupling, is 
expected to cause only mild experimental difficulties. 

As an example of Berry's phase for half-integer spins 
we have chosen 5 — 3/2. Our predictions can be ex- 
perimentally verified on ^°^IIg ground state atoms. For 
positive values of the quadrupole to dipole coupling ratio 
A ~ 2, the state 5 = 3/2, m = —1/2 exhibits a large 
Berry's phase in strong contrast with the 5 = 1/2, m = 
— 1/2 state behaviour. As seen in Sec.V and Appendix 
B, taking advantage of this 5-dependence enables one to 
produce maximum entanglement between three initially 
non-correlated 1/2-Qbits in the state m = —1/2. 

As Sec. VI demonstrates, our goal aiming at the explo- 
ration of AA cycles by realizing parallel transport of spins 
5=1 looks achievable. Thanks to precisely adjusted 
time dependences of the light-shift Hamiltonian parame- 
ters, (intensity, helicity of the beam, polarization rotation 
speed and magnitude of the magnetic field colinear to the 
light beam), it is possible to satisfy the parallel transport 
condition at any instant t of the cycle. As a result the 
trajectory followed by the spin-1 quantum state during 
its cyclic evolution is a closed loop composed of geodesic 
segments drawn on the spin-1 state space, i.e. the com- 
plex projective plane CP^. This provides new possibili- 
ties for exploring empirically the non-trivial geometrical 
properties of this four-dimensional space. 

From a physical point of view, the main interest of 
the AA geometrical phase is that it could be generated 
by using as a 5 = 1 candidate the F = 1 hf level of the 
ground state of ® ''Rb, which is the building block of one of 
the most popular "optical" crystal. As for Berry's phase 
there is a large flexibility available on the cycle duration 
provided correct scaling is applied to the other parame- 
ters, beam intensity and magnetic field. The advantage 
of the spin-1 AA phase compared to Berry's is twofold. 

• a. All limitations coming from the necessity of get- 
ting rid of the non adiabatic corrections disappear. 

• b. The instability problem coming from the use 
of effective B, E fields are much easier to control. 
This follows from the remarkable fact that the atom 
loss during the cycle duration Tc depends only on 
the value of the AA geometrical phase and not upon 
Tc (see equation ([551) 1. For a typical value of /3aa 
the percentage of atom loss per cycle is about 15%. 
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In the present paper we have dealt, for sake of simplicity, 
only with cycles involving a and A as time-dependent 
coordinates. The construction of a parallel transport 
Hamiltonian can be extended to an arbitrary closed cir- 
cuit of CP^, using the results of reference |28| . 



Appendix A: Entanglement of three non-correlated 
one-half spins using Berry's quantum cycles 

Adapting the method of reference [i| , we introduce the 
three non-corellated spin states with M = X^i'^i '^li = ^: 

(Al) 

(E"*^^^ and <I>(^^ being obtained by a circular permutation 
of the jTLi. These three states form an orthogonal basis 
for the set of three one half-spin states with M = ^ . The 
next step is to construct three orthogonal eigenstates of 
with S = X]i=i which are linear combination of the 
Using the rules of adddition of quantum angular 
momenta, one finds that the possible eigenvalues of S^, 
h'^S^S + 1), correspond to 5 = | and S = ^. There is 
a unique way to construct the state S — ^: one applies 
the operator S- = Sx — iSy upon the state S = M = |, 
i.e. ^'3 3 = |i) ® |i) (g) One obtains immediately: 

^3 1 ^ ^($(1) +$(2) +$(3)). (A2) 

This state is invariant under all permuations of the 3- 
spin states. One has now to construct two orthogonal 
states with S = M = i, denoted ^I'" i , (n— 1,2), which 

2 2 

differ by their symmetry under the permutations of the 
three spins. Ignoring for a moment the orthogonality 
condition, it is easy to obtain two such states, linearly 
independent, $(1-^) = - $(j))/v^ with j = 2,3. 

By introducing in ^^^J) the explicit expression of the 
states ^(■'^ one finds, by applying the rising operator 
S+ = Sx+ iSy, that indeed, S+^^i'-J) = 0. It is then 
easily seen that two orthogonal states vj/^ i are given, up 

2 2 

to a normalization factor, by the sum and the difference 

$(1,2) ^ $(1,3). 



^1 ^ ^ J_(2$(l)-$(2)_$(3);^2 ^ ^ J_($(2)_<J,(3)-)^ 
2 2^ '22^ 

(A3) 

It is now a matter of simple algebra to write the non- 
correlated state as a linear combination of the three 
above angular momentum eigenstates: 

$(1) = ^(5-3 1 V2*\ 1) (A4) 
V3 ^= ^5 

An important feature of the three states 

^'3 1 , ^'1 1 , ^1 j_ is their symmetry properties un- 

22 22 22 

der permutations of the three spins. Let us introduce 



the fully symmetric operator S = Si=i Pi where 
Pi is one of the 6 possible permutations of the three 

TOi. It is then easily verified that: 5 4*3 1 = vps 1 

2222 

while S "^^l 1=0. On the other hand, the permuta- 

2 2 

tion (23) applied upon the two S — ^ states gives: 

(23) *i 1 = ^-i 1 and (23) ^'i 1 = 1 ■ 

22 22 22 22 

We are going to study the adiabatic evolution of 
the three non-correlated, 1/2 spin states governed by 
Hsit) = H{B{t),'E{t)), which looks formally like the 
quadratic spin Hamiltonian of Eq.([T]) discussed exten- 
sively in the present paper, but with a crucial difference 
lying in the fact that S is meant to be the total spin op- 
erator S = X]i=i ^i- The Hamiltonian H3{t) is invariant 
under all the the permutations of the three spins and, as a 
consequence, all its non-diagonal matrix elements taken 
between any pair of the three above states are vanish- 
ing. The above result can be extended to the set of three 
states ^'i M> *i M' ^^^^ 1^1 = l^^'l = 5 obtained by ap- 
plication of the raising (lowering) operators S± — Sx^Sy, 
also permutation invariant. We are then lead to the con- 
clusion 1] that the three states i m^^\ m' behave vis 
a vis the Hamiltonian H^it), as if they were associated 
with isolated spins S. 

Our Berry's cycle would be organized in a way similar 
to the one we described in section IV of the present paper, 
but with one difference: the precession Euler angle (p has 
to be among the cyclic parameter in order to have a non 
vanishing Berry's phase for 5 = i. Otherwise it would be 
impossible to achieve a maximum entanglement. As an 
illustration, we shall consider situations where during the 
(/j-cycles 6 and A have predefined fixed values obtained by 
adiabatic ramping processes, analogue to those discussed 
in [IJ. It is then convenient to introduce the difference 
A/3{X,9) between the 5* = | and S = ^ Berry's phases: 

A/3(A,0) =cos0 (p(|,i;A)-i)2rv^. 

At the end of the Berry's cycle, the initial non- 
correlated state <E>*^^-', written as a linear combination 
of angular states behaving as isolated spin systems, has 
evolved into the following state: 

$gl,(A) = exp(z/3(|i;A)-i=*3 1 -I- 

exp(z/?(i,i;A)vI'i 1. (A5) 
2 2 22 

We rewrite the right-hand side of the above equation in 
terms of the non-correlated states and factor out 

the overall phase x = l^i^,^,^) in order to exhibit the 
Berry's phase difference A/3 = /3(|, i, A) - /3(i, i). We 
get the final expression for our candidate for a three one 
half-spin entangled state: 

^^bpW = ^^((l + 2exp(-iA/3))<l>(l) + 

(1 - exp(-iA;3))($(2) - $(3)). (A6) 

To achieve a maximum entanglement one must impose 
the equality of the absolute values of the two mixing co- 
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efficients, |1 + 2exp(-iA/3)| = |1 - exp(-iA^)|. It is 
easily found that this implies the condition A/3 — 
To get explicit results, we have chosen the typical case 
riip = 2 and 9 = Using equation (jB2p of Appendix 
B, one gets the following value for X,nax — —0.784562. 
This negative value is welcome because it allows one to 
avoid the near-crossing of the two levels f (|, |, A) and 
f(|,5,A) which occurs for A > 1. For = 2 and 
9 — ^ the S — ^ Berry's phase takes the simple value: 
/3(^,^) = 7rmod(2 7r), since we know A/3, this leads to 



) = — -I mod(2 tt) and to the final form 
<&L'kA™a.) = -^(-exp(i^)|-i)®|i)®|i) 
+ expH^)|i)$5|-i)(»|i) 
-expH^)|i)®|i)®|-i)).(A7) 

The initial non-correlated state has been trans- 
formed, at the end of this specially designed Berry's cycle, 
into a correlated state with a maximum entanglement. 



reference [29| in order to make them compatible with 
the notations of the present paper. We begin by the re- 
duced Hamiltonians Heven and 'Hodd connecting states 
with (—1)3/2-™ ^ _|_-[^ respectively: 



3A 



-He 



V 2 4 2 



V3\ 
2 

_A _ L 
4 2 



7A 



(Bl) 



From them one gets readily the explicit values of £{m, A) 
I and m — ^. 



and p(m. A) for m 

£(3/2, A) - 
p(3/2,A) = 



X 



1/2 



2VA2 - A + 1 
f(l/2,A) = i (4VA2 + A + 1 + 
A + 2 



5A - 2 



p(l/2,A) 



2VA2 - A + 1 



1/2 



(B2) 



Appendix B: Berry's phase quantum cycles 
generated by S= 3/2 non linear Hamiltonians 



The formulas for m < are easily obtained from the 
reflexion law derived in ll section III: 



We give here basic formulas for the Berry phases rel- 
ative to a spin S — ^. They have been adapted from 



£{—m, A) — ~£{m, —A) , p{—m, A) — —p{m, —A). 
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